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I. Introduction 



In previous work |7j the authors considered the asymmetric simple exclusion pro- 
cess (ASEP) where particles are confined to the nonnegative integers Z^ = {0, 1,2,.. .}. 
Each particle waits exponential time, and then with probability p it moves one step to 
the right if the site is unoccupied, otherwise it does not move; and with probability 
^ ■ g = 1 — p a particle not at moves one step to the left if the site is unoccupied, 

Q . otherwise it does not move. For n-particle ASEP a possible cofigurations is 

00 : 

O. X = {xi,...,x„}, (0 < Xi < ■■ ■ < x„). 

^' 

O ■ The Xi are the occupied sites. We denote by Xn the set of possible configurations 

for n-particle ASEP, and by pn(x, y;t) the probability that at time t the system is 
in configuration x given that initially it was in configuration y. (We shall drop the 
subscript 



H ' In [7] a formula was found for this probability. It was the sum over the Weyl 

■ ■ group B„ of multiple integrals. (For ASEP on Z it was a sum over the permutation 

group §„ [SIE].) 

Here we consider the ASEP on Z+ with an open boundary at zero. The stationary 
measure for ASEP on the finite lattice [1, L] or on the semi- infinite lattice Z"*" with 
boundaries connected to reservoirs has been the subject of much research starting 
with Derrida et al. [1\. We refer the reader to the recent work of Sasamoto and 
Williams [1| for an up-to-date account of these developments. Here we consider the 
time-dependent properties of ASEP on Z"*" with an open boundary. Specifically, the 
point is connected to a reservoir where a particle is injected into site from the 



reservoir at a rate a, assuming that the site is empty, and a particle at site is 
ejected into the reservoir at a rate /3. Now the number of particles is not conserved 
and for ASEP with open boundary the configuration x may lie in Xn while y may lie 
in Xm with m y^ n. 

We find an infinite tri-diagonal matrix with operator entries in which the Laplace 
transforms of the probabilities can be read off from the entries of the inverse matrix. 
When either a = or /3 = the matrix is triangular and so the inverse can be 
computed more explicitly. The result is obtained by solving a system of differential 
equations for the probabilities. The final formulas involve inverses of operators with 
kernels the Laplace transforms of certain p(x, y;t) obtained in [7]. 

There are two special cases in which the results are more explicit. For TASEP 
with p = 1, the inverse operator is computable in terms of p(x, y;t) itself, and the 
probabilities are given in terms of certain determinants. For SSEP {p = q) and general 
a and /3 we find formulas analogous to the ones described above for the probability 
that sites Xi, . . . , a;„ are occupied. This is for infinite systems as well as finite ones. 

We state the formulas for p(x, y; t) in the Appendix. 

II. Statement of results 

We denote by £"„ the Lebesgue space L^{Xn)- From the fact that 

for each y, it follows that the operator on £„ with this kernel is bounded with norm 
one. We denote the Laplace transform of p(x, y; t) by p(x, y; s)\ 



POO 

P(x,y;s) = / p(x,y;t)e" 
Jo 



St 



dt. 



The operator with this kernel is bounded on Sn with norm at most (Res)^^. We 
denote it by Ln(s)|j In the results stated below it is tacitly assumed that Res is 
sufficiently large. 

Now for ASEP with open boundary at zero, we define Pn(x; t) to be the probability 
that the system is in configuration x G A'„ at time t. (We shall usually drop the "x" 
in the notation, and do not specify an initial configuration.) We denote its Laplace 
transform by Pn{s). 

We define vector functions 

Pis) = (P„(s))„>o, P(0) = (P„(0))„>o, 



^When n = we interpret p(x, y;i) as 1, and so Lo{s) is multiplication by s ^. 



belonging to the direct sum Xl^o^"- 

We define operators An : Sn-i — > Sn and Bn : Sn+i — ^ Sn by 

{AnF){xi,...,Xn) = S{Xi)F{x2,...,Xn), (2) 

(5„F)(a;i, . . . , x„) = (1 - 5(xi)) F(0, xi, . . . , x„),| (3) 

Then we define matrices 6, L{s), A, B, with operator entries, acting on ^^q^"- 
The first is diagonal with (ra, 'n.)-entry multiplication by S{xi), the second is diagonal 
with n, n-entry L„(s), the third is subdiagonal (one diagonal below the main diagonal) 
with n,n — 1-entry An, and the last is superdiagonal (one diagonal above the main 
diagonal) with n,n + 1-entry Bn- 

Theorem 1. With this notation we have 

P{s-a)= (l-L{s){{a-(3)S + aA + l3B)) L{s) P{0). (4) 

There are expressions for the entries of the inverse operator as infinite series of 
products. But when either /3 = or a = the operator has only one subdiagonal 
or one superdiagonal and each entry of the inverse is a single product. We state the 
results as recursion formulas. Define 

Mn{s) = {I-{a-/3)Ln{s)6)-\ (5) 

Corollary 1.1. Suppose /3 = and that initially there are k particles at y G A'fc. 
Then 

Pkis-a) = Mkis)Lkis)6y3 

and when n > k 

Pn{s -a) =a M„(s) L„(s) A„ P„_i(s - a). 

Corollary 1.2. Suppose a = and that initially there are k particles at y G Xk- 
Then 

Pkis) = Mki-s) Lkis) 6y, 

and when n < k 

Pni-S) = P Mni-S) Lni-S) BnPn+lis). 

In connection with the corollaries we show the following. 



^When n = we interpret S{xi) as zero. In particular Aq — and Bq F = F{Q). 
■^This is interpreted as s~^ when k = 0. 



Remark 1.1. The operators appearing in the inverses can be replaced by lower- 
dimensional ones. This will be useful for computation. Define 

A-^ = {{xi, . . . ,a;J e A'„ : xi > 0}, C = L^X^^), 

and then operators: 

Ll-M : C-i ^ C-i with kernel p((0,x), (0,y); s), 
Ll,n-ii.s) ■ ^n-i -^ £n with kernel p(x, (0,y); s), 
^n-i,n(s) : £n ^ £n-i with kernel p((0,x),y; s). 

(a) The operator Mn{s) Ln{s) An : £n-i — >■ £n in Corollary 1.1 is equal to 

Lln-i{s) {I - {a - (5) Ll_,{s))-' RnS 

where Rn-i '■ £n-i — ^ £n-i is the restriction operator. 

(b) The operator Mk{s) Lk{s) : Sk -^ Sk in Corollaries 1.1 and 1.2 is equal to 

Lkis) + (a - /3) L'is) (I -{a- /3) LUi-s))-' LUJs). 



Remark 1.2. In the special case of TASEP when p = 1 we have the simplification 

{I -a LO (s))-i = I + a Ll{s - a). 

In the case of SSEP {p = q), even for infinitly many particles, there are formulas 
for correlations that are no more complicated when both a and (3 are nonzero. For 
X = {xi, . . . ,Xn} € Xn we define \E'n(x; t) to be the probability that sites xi, . . . , x„ 
are occupied at time t. We denote its Laplace transform by \E'n(-s) = ^E^nl^; s) and 
introduce the vector functions 

$(s) = ($„(s))„>o, ^(0) = (^„(0))„>ol 

Let the operators Ln{s) and An, and the matrices S, L{s), A with operator entries 
be the same as beforecl 

Theorem 2. We have, 

$(s) = (/ + L{s) {{a + (3)6-a A)) L{s) ^(0). 



''The / here is the identity operator on S^^i while the / in ([5|) is the identity operator on f„. 

^We define ^o(t) = 1, and so $o(s) = s"^ 

^Now we define £„ = L°°(A'„), and observe that by H]) and the symmetry of the kernel we have 
J2y£X P(^7y;0 = 1 fo'" each x e A",;. It follows that L„(s) is a bounded operator on this f„ with 
norm at most (Res)"'^. 



We now set 

M„(s) = (J+(a + /3)L„(s)5)-\ 

Corollary 2.1. For n > 0, 

§„(s) = a Mn{s)Ln{s)An ^ n-l{s) + Mn{s)Ln{s) ^„(0). 

Corollary 2.2. In the case of Bernoulli initial condition with density p, 

§„(s) = a Mn{s)Ln{s)An ^n-i(s) + s"^ p" M„(s) 1, 

where "1" is the constant function on Xn. 
Corollary 2.3. When initially no sites are occupied, 

§„(s) = aM„(s)L„(s)A,$„_i(s). 

The analogue of Remark 1.1 holds here. 

From the abstract formulas, Theorems 1 and 2 and their corollaries, we derive 
some concrete results. 

Suppose, in ASEP, that at time zero there is a single particle at y. From Corol- 
laries 1.1 and 1.2 we show: 

q: = 0: When p > q, with probability 

^ (3{q-^-i)-y 

p- q + l3 

the particle is never ejected. When p < q, with probability one the particle will 
eventually be ejected. The expected value of the time this occurs is infinite when 
p = q, and when p < g it is 

y + q/P 

q-p 

/3 = 0: With probability one a second particle will eventually be injected. The 
expected value of the time at which this occurs is 

1 ^ e+(«)-^ 



a q{^+{a)-l)-a 

where 

e+(a) = !(« + ! + y/(a + iy-Apq) . 



Combining Remark 1.2 with the determinant formula (!24|) for p(x, y; t) in TASEP 
it is practical to compute some exact results. Denote by 'P„(t) the probability that 
starting with no particles at time there are exactly n particles at time t. Then 



Vo{t) = e 



+ 



Viit) 



^2(t) 



a 
2(1 -a] 
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1 — a 
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These are for a ^ 1. When a = 1 there is no singularity; we take a limit and the 
formulas simplify. 




Figure 1: Plotted are the TASEP probabilities 'P„(t) for n = 0,1,2,3 with a = 1. 
Increasing n moves the maximum to the right. 



For SSEP we obtain the following consequence of Corollary 2.2. Define 

ANit)=J2ivAt)-vM). 



x>0 



the net number of particles that have entered the system at time t. (Which may 
be negative.) With Bernoulh initial condition when {f]x{0)) = p, we show for the 
expected value that 



(AiV(t)) 



V TT a + /3 



By a laborious computation (not included) we can show that the second moment of 
AiV(t) is finite. 

In the special case where initially there are no occupied sites, we use Corollary 2.3 
and present a not completely rigorous (to say the least) argument that 

(AiV(t)^)^- " t ast->oo. 
TT {a + f3)'^ 

Combining this with the first moment asymptotics when p = we conclude that the 
variance of AN{t)/t^/'^ tends to zero as t — > oo|3 The derivation is quite long, but 
the result with the precise constants came out so nicely in the end that we could not 
resist including it. 

III. Proofs of the results 

A. ASEP 

Proof of Theorem 1: The probability p(x, y; t) for n-particle ASEP on Z+ is the 
solution of the differential equation 

, N 

— u{X;t) = ^ lpM(xj- 1) (1 -5{xi -Xi-i - 1)) + qu{xi + l) (1 - 5{xi+i ~ Xi - 1)) 
j=i 

-pu{xi) {l-8{xi+i-Xi-l))-qu{xi) {l-5{xi-Xi^i-l)) +[qu{xi)-pu{xi-l)] 5{xi) 
that satisfies the initial condition 

u(x; 0) = 5y(x). 

(In the ith summand entry i is displayed and entry j is Xj when j ^ i.) 

If we denote by Q„ the operator given by the right side of the equation, then 
p(x, y; t) is the kernel of e*®". Thus the equation is 

dpr 



dt 



Q„p„(x; t). 



^This result was obtained in [5] in the case /3 = and p = 0. 
^As predicted in l2j in the case /3 = 0. 



where we have written pn(x; t) for p(x, y; t). 

For open ASEP we write the corresponding probabihty as P„(x; t) (and do not 
specify any initial condition). The equation for P„ is 

dPn 

— -T— l^Xi, . . . , Xfi . Zj tglfj J jil^Xi, . . . , Xji . Z) 

+a 6{xi) Pn-i{x2, ...,Xn: t) - a (1 - 6{xi)) Pn{Xi, . . . , x„ : t) 
-/3 6{xi) Pn{xi, ...,Xn: t) + /3(1 - 6{xi)) P„+l(0, Xi, . . . , x„ : t). 

Define $n(x; t) = e"* P„(x; t). In terms of the operators An and i?„, defined in ([2]) 
and (IHl), and 5 = S{xi), the equation for $„ becomes 

d^ (t) 

— ^ = Qn $„ + (a - /3) 5 $„ + a A„ $„_i + /3 P„ $„+i. (6) 

The equation and initial condition are satisfied if 

$n(t)= / e(*-")2n(^(a_/3)5$„(M) + aA„<l>„_i(n) + /3P„<l>,,+i(n))dn + e*2"$„(0). 

We use the fact that the Laplace transform of a convolution is the product of the 
Laplace transforms. Recall that the kernel of Ln{s) is p(x, y; s), which is the Laplace 
transform of the kernel of e*^" . In other words, the Laplace transform of the operator 
e*®" is Ln{s). So taking Laplace transforms in the last equation gives 

$„(s) = L„(s) (^(a-/3)<5$„(s) + aA„$„_i(s)+/3P„$„+i(s)) +L„(s)<l'„(0). 

If we now introduce the vector functions 

$(s) = ($„(s))„>o, $(0) = (<l'„(0))„,>o, 

and the operator matrices 5, L{s), A, B defined earlier we see that the system may 
be written as 

$(s) =L{s) ((a-P)5 + aA + PB\^s) + L{s)^{'d). 
Since $(s) = P{s — a), this gives the statement of Theorem 1. 
Corollaries 1.1 and 1.2: We write the operator inverse in (j4]) as 
'l -{M{s)L{s){aA + pB))] M{s)L{s), 



where M{s) is the diagonal matrix with operator entries M„(s). When /3 = this 
equals 

{I-aM{s)L{s)A)'^M{s)L{s). (7) 

The operator matrix M{s) L{s) A consists of one subdiagonal, with m,m — 1-entry 
Mm{s) Lm{s) Am- Therefore the n, n'-entry of the inverse {n' < n) is 

a"-"' Mn{s)Ln{s)An ■ ■ ■ M„/+i(s)L„,+i(s) A„/+i, 
where for n = n' this equals /. So the n, n'-entry of (171) is 

a"~"'M„(s)L„(s)A„---M„/+i(s)L„/+i(s)A„'+iM„/(s)L„/(s). 
Thus, 
Pnis -a) = J2 «"""' Mnis)Lnis)An ■ ■ ■ M„,+i (s)L„,+i(s) A„,+i Mn'i-s) Ln'{s) P„.(0), 

n'<n 

(8) 

When a = the matrix M{s) L{s) B consists of one superdiagonal, and we obtain 
similarly 

P„(s) = ^ /3"'-" M„(s)L„(s)5„ ■ ■ ■ M„,_i(s)L„,_i(s)5„,_i M„,(s) L„,(s) P„,(0). 

n'>n 

If initially there are k particles at y G A/c then in both cases Pn'(O) is nonzero 
only for n' = k. The formulas become 



Pni-s - a) = a" M„(s)L„(s)A„ ■ ■ ■ Mk+iis)Lk+iis)Ak+i Mk{s) Lk{s) 5, 



yi 



Pnis) = /3'=-" M^{s)L^{s)B^ ■ ■ ■ Mk-i{s)Lk-i{s)Bk-i M^i-s) Lk{s) 6y, 
and the corollaries follow. 

Remark 1.1: For (a) we use the fact that because An has the factor S, 

Mn{s)Ln{s)An = {I - a L„(s) 5)-^Ln{s) SAn = L„(s) 5{I-a5 Ln{s) 5)-^ An. 

Since (A„/)(xi,x) = 5(xi) (-R„_i/)(x) and x G £n-ii we obtain statement (a) in 
different notation. For (b) we use 

Mfc(s) = {I-{a-l3) Lk{s)6)-' = I + {a - /3) Lk{s)6 (I - (a - (3) Lk{s)6y' 

= /+(«-/?) Lkis)6 (I -{a- /3) 5Lkis)6)-\ (9) 



Thus, 

M,(s)Lfc(s) = Lk{s) + {a-/3) Lk{s)6 {I - {a - /3) 6Lk{s)6)-'6Lk{s), 
and statement (b) follows. 

Initially a single particle: Let Vi{y;t) = J^xyo^^i^^y^^) denote the probability 
that, starting with one particle at y, we still have one particle at time t. Denote its 
Laplace transform by Vi{y; s). 

We begin with the case /3 = 0, so Vi{y; t) is the probability that no new particle 
has been injected by time t. 

By Corollary 1.1 the Laplace transform Pi{x, y; s—a) is equal to {Mi{s)Li{s)Sy){x). 
By Remark 1.1(b) this equals 

p{x, y- s) + ap(a;, 0; s) (1 - ap(0, 0; s))"^ p(0, y- s). 
Then using X]x>o P('''' V'^ ^) ~ ^~^ *^ obtain 



Vi{y]s-a) 



1 + 



ap(0,z/;s) 



l-ap(0,0;s) 
From formula ( 123|) for p(x, y\ t) in the case n = 1 we compute that 

1 i+{s)-y 



P(0,y;s) 



qUs) 



(10) 



where 



^+{s) = 1 (s + 1 + ^{s + lf-Apq)- 



this is the solution of e{^) = s with positive square root when s > 0. Thus 



Vi{y]s-a) 



1 + 



ai+{s)-y 



q{i+{s) 



a 



'IV. 



The denominator in the brackets is nonzero for s = a (and positive for s > a), from 
which we conclude that 

V^{y;t)dt = V,{y;0) = -+ ,, ^"^ \^ . 

a ?(^+(a)-l)-a 

It follows that with probability one a second particle will eventually be injected, since 
V{y; t) — > as t — i- oo, and the integral is the expected time at which it occurs^ 



^If T{y) denotes the time when a second particle is injected, then Prob (r(y) > t) = Vi{y\t), 
from which the statement follows. 



10 



This was for /3 = 0. For a = 0, Vi{y; t) is the probabihty that the particle has not 
been ejected by time t. We use Corollary 1.2 (and Remark 1.1), and formula (ITT]) is 
replaced by 



Vi{y;s) 



We compute that as s — )■ 0, 



(3 Us) 



'Piiy^s) ~ - 



g(e+(s)-l)+/3 
/3(g-i-l)-^ 



p- q + (3 



if p > q, 



Vi{y;s) 



{2y + l/f3) iip = q, 



^1 y; s) -> if p < g. 

q-p 

Applying the Tauberian theorem we deduce 

R(n-^ _ i)-y 

]imr,{y;t) = l~'^ -j— iip>q, 

t^x> p — q + p 



Vi{y;t) 



1 



2y + 



1 



t ^'"^ as t — > oo if p = g, 



/2^L '^ /3J 

V,{y-t)dt = ^^^^^ iip<q. 
q-p 

When p > q the limit on the first line is the probability that the particle is never 
ejected, li p < q, then with probability 1 the particle will eventually be ejected. The 
expected value of the time at which it is ejected is infinite when p = g, by the second 
line, and is given by the integral on the last line when p < q. 

B. TASEP 



Remark 1.2: To compute (/ — a L^{s)) ^ we use the fact that for k > 1 the kernel 
of LO(s)^ is 

J2 P((0, x), (0, zi); s) p((0, zi), (0, Z2); s) ■ ■ ■ p((0, z,_i), (0, y); s). 

zi,...,Zk_i£XZ 

The summand is the Laplace transform of the {k — l)-fold convolution 



p((0, x), (0, zi); ui) p((0, zi), (0, Z2); U2) ■ ■ ■ p((0, z„_i), (0, y); Uk^i) du. 



Ul-\ |-Mfc_l=t 



11 



When p = 1, if the left-most particle begins at and ends at then it was always at 0. 
It follows that these probabilities have the semigroup property. Thus after summing 
over zi, . . . , Zfc_i the integral becomes 

/ p((0,x),(0,y),t)rfu=-^-— p((0,x),(0,y),t). 

Jui-\ — i-«fc_i=t y'^ -*-/• 



Therefore the kernel of of L°(s) is 



lis)' ie 

OO J-fc— 1 



/'OO 4.K— I 

It follows that the kernel of ^^i(a-^°(s))'^ is 

a / e-^*e"*p((0,x),(0,y),t)rft = «p((0,x),(0,y);s-«). 



This gives (J — a -^° (s)) ^ = 1 + a L^{s — a). 

Once we have this result (and Corollary 1.2, Remark 1.1, and the determinant 
formula ^2M\ for p(x, y;t)), we have all the ingredients necessary to compute Vnif) 
for small n. This is the probability that starting with no particles at time there 
are exactly n particles at time t, which equals XlxeA" 'Pn{'^',t). The computation of 
products of kernels involves summing geometric series. 

C. SSEP 

Proof of Theorem 2: A state of the system is a function rj : Z+ — )■ {0, 1} where 
rj^ = 1 means site x is occupied and ri;^. = means site x is not occupied. Recall that 
we defined \l/„(xi, X2, . . . , x„; t) as the probability that sites xi, . . . ,Xn are occupied 
at time t. Thus, 

^n{xuX2,...,Xn; t) = (r^^i (t) ■ ■ ■ r/^„ (t)) . 

The Markov generator C of ASEP on Z+ with an open boundary at zero is given 
mby 

Cfiv) = a (1 - vo) UXv') - fiv)) +PVo ifiv') - f{v)) 

OO 

+ J2[pVk{l- %+i) + g (1 - V,) %+i] {fiv'''^') - fiv)) . (12) 

Here / is an M-valued function that depends on only finitely many sites, and 

r 1 71 if T - A- f "^^+1 ii x = k 

l^ ^" II a; ^K, 1^ ^^ if x^A;,A; + l. 

12 



For SSEP the sum in (fT2|) equals 1/2 times 

oo 

J^lvk (1 - vk^i) + (1 - Vk) vk+i] UXv"'"^') - fiv))- 



k=0 



When rjk = rjk+i the first factor equals 0; otherwise it equals 1. Since the second 
factor is zero when rj^ = %+!, we can ignore the first factor, and we get 



x:(/(^'''^')-/(^))- 



k=0 



For the correlations we are interested in f{ri) = rj^^ ■ ■ ■fjxny so the kth summand 
equals zero unless either k = Xi for some i ot k = Xi — 1 for some i. (We'll see that 
these cannot both happen for a nonzero summand.) 

Suppose first that k = Xi. If Xj+i = Xj + 1 then k + 1 = Xj+i and the substitution 
r] — )■ r/'^''^+i applied to /(?7) just interchanges t]^^ and fjx.+i- Therefore the kth summand 
is zero. It follows that for a nonzero summand we must have Xi+i > Xi + 1 = k + 1, 
and the substitution rj — )■ rj^^^^^ only affects the ith factor in f{r]). Therefore the 
summand equals 

(r]x^ ■ ■■Vx,.^ Vx,+iVx,+i ■ --Vx^-Vxi ■ ■■Vx,.^ Vx, Vx,+^ ■ ■■Vx„) (l-(5(xi+i-Xi-l)). (13) 
Similarly if A; = x, — 1 the summand equals 
(^^■1 ■ ■ ■ ^x,_i Vx,-i ^x,+i ■ ■ ■ r]x„-Vx^ ■ ■ ■ Vx,_, Vx,, Vx,+i ■ ■ ■ Vx„) {l-6{xi-Xi_i-l)). (14) 

(This is to be multiplied by 1 — 6{xi) when i = 1.) 

li k = Xi for (TT^ and k = Xi' — 1 for f lT^ then i' = i + 1, so Xj+i = Xi' = Xi + 1 
and (IT^ zero, and Xj/ = Xj + 1 = Xj/_i + 1 so flT^ is zero. Thus for the fcth summand 
to be nonzero either k = Xi for some i ot k = Xi — 1 for some i, but not both. 

It follows that for SSEP the expected value of the sum in ( 1T^ is equal to 1/2 
times the sum over i of the expected values of the sum of ( ITSjl and ([H]) . This equals 

Qn{{Vxi---Vxn))- 

Adding what we get from the first two terms of (II 2p we find that the differential 
equation for \E'„ = \E'n(t) is 

— *„(Xi, X2, . . . , X„) = Qn ^„(Xi, X2, . . . , X„) 

at 
-{a + (3) ^„(xi, , X2, . . . , x„) 6{xi) + a ^„_i(x2, . . . , x„) 5(xi)0 



^"Observe that the constant functions ^„ = (a/(a + /3))" satisfy the equations. Thus the Bernoulli 
measure with density p = a/{a + /3) is stationary. This also follows from the results of [4] in which 
stationary measures were determined for general ASEP on Z+. 
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Except for the change a — /3 -^ — (« + /3), this is (|6]) without the Bn terms. Therefore 
to complete the proof of Theorem 2 we need only make this change as we go through 
the rest of the proof of Theorem 1, which we need not do. 

Corollaries 2.1 and 2.2: Just as ([8]) is obtained from Theorem 1, we obtain now 
for n > 0, 

n 
^n(s) = J2 """"^ Mrr{s)Lr^{s)An ■ ■ ■ M,.+i (s)Lfc+i(s) Afc+i Mfc(s)Lfc(s) ^fc(O). 
fc=0 

Corollary 2.1 follows. For Corollary 2.2 we have \E'n(0) = p", and we use the fact that 
L„(l) = XlyeATn P(^' yj ^) ~ ^^^- Corollary 2.3 is the case p = of Corollary 2.2. 

A.N{t), the net number of particles that entered the system: We assume 
that we have SSEP with Bernoulli initial condition. Corollary 2.2 when n = 1 gives, 
for the Laplace transform of (//^.(t)), 

Sr)(s) = $i(x; s) = s-^ (pMi(s) 1 + aMi(s)Li(s)Ai 1). 
By Remark 1.1(a) applied here, we hava^^l 

Mi{s)L,is)A,l = -^^^^^l^ (15) 

l + 7p(0,0;s) 

where we set 

7 = a + /3. 

Similarly, from i^ we obtain 

M,(.)l = l-7: ^^"'°'^^ 



l + 7p(0,0;s) 
Combining the two gives 



s 



.^ l + 7P(0,0;s)^^ ' 



Subtracting p/s from both sides and summing over x > we get for the Laplace 
transform of {AN{t)) 



s2 l + 7p(0,0;s)' 



^^ Since the definition of Mk(s) is different now we must replace a — /S by — (a + /3) when using 
the remark. 

14 



From ( ITOj) we have for SSEP 

2 [2 

p(0, 0: s) = , ~ W - as s -> 0. 

Hence 

(AiV(t))(g)^ ^-7P ^-3/2 ass^O. 

By the Tauberian theeorem this imphes 

2 a-7P^i/2 



(AAr(t))~ W-:: — ^t^^ ast^oo. 

TT 7 



As stated in the last section the second moment {AN{t)^) is finite. To show this 
we use that the second moment is equal to 

lim Yl ((^-i(^) - ^-i(0)) (^-2(^) - V.M))- 

xi,X2<N 

We can show that the sum is a polynomial of degree two in p, and that each of the 
three coefficients of the powers of p has a limit as A^ — > oo. The argument is quite 
involved, and we do not include it. 

The second moment of A.N{t) when p = 0: Now it is mainly a question of 
determining the asymptotics of '^xi<x2(''l^i^^2) ■ The Laplace transform of \E'2(x; t) = 
(Vxi VX2) is given by Corollary 2.2 as 

$2(x; s) = a" S-' M2{s)L2{s)A2Mi{s)Li{s)A^ 1. 
We computed Mi{s)Li{s)Ai 1 in flT^ . Combining this with Remark 1.1(a) gives 



^2{^-s) = a's-'LUs){I + lL%s))-'- 



Po 



l + 7Po(0)' 

where po is the function x — t- p(a;, 0; s). 

To obtain the Laplace transform ^^ ^^ ('7x1^x2) we sum over x e ^"2. If we recall 
that the kernel of Lg ^(s) is p(x, (0,y)) and that the sum of this over x G ^"2 is s~^ 
we see that the desired sum is the inner product of the remaining function with the 
constant function s~^. Thus, 

E (^fa^) = , f%].A {i+iLUs)r%A). (16) 

x77^2 i + 7Po(o) 
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(Here we use the fact that po belongs to L^ and that L^{s) is a bounded operator on 
this space.) 

What follows is not rigorous. We want to rescale (/ + 7L°(s))~^ as s — )■ 0, and 
we refer to formula fl23|) given in the Appendix for p(x, y;t). After taking Laplace 
transforms in the case n = 2 we find that L^{s) has kernel 



,„, (2") Vc. ^c. "»-£&) -£(6) 6& 

where Cr is a circle with radius R with R large. (Some integrals are taken over two 
pairs of different contours and the results averaged.) To begin with, s is so large that 
taking Laplace transforms under the integral sign is valid. 

If we ignore the poles of the A^j we can move both contours to the unit circle Ci. 
Then the range of e{^i) + ^(^2) is [—4, 0], so we may take any s > 0. As s — > the 
main contribution comes from a neighborhood of ^1 = ^2 = 1 because the denominator 
vanishes there when s = 0. If we set ^1 = e™^, ^2 = 6**^^ then the integral with its 
factor becomes to first order 

I r r g«(«<T(2)^-i'2y) 






^ „ , , A^(e-V2..i .V2..2) - dVidV2. 



This becomes, after the scaling x — >■ a;/v2s, y — ?• Z//v2s 

^2-, 



^, , A,(e*^2^''\e'^2s.2)r ^dvrdv2, 

which acts on functions on M^. 

Each A(j has absolute value 1 on Ci x Ci, and each A^j^e''^^'"^ , e^^'"'^) has limit 1 
as s — !■ except when vi = V2 (mod 27r). Thus we replace the above by the approxi- 
mation 

— ^ , — / / -; 5 7^dvidv2. 

271^^2^ kU l + vl + vl 



^2 

This depends only on cr(2). If we use the fact that the denominator is even in each 
Vi, and that V-i = —Vi, we see that the sum over a G B2 of the integrals equals 

1 f f e*^2 '•^"^^ + e*"^ ^^"^^^ 2 f f e^'^^i^~^2y) 

dvi dv2 -\ ^_- / / :;— — 5— — 2 '^'^1 '^^2- 



7i'^V2sJrJr 1 + vj + vi 7^^^^/2^ MJul + vj + u 
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In the first integral we integrate first with respect to fi, and we obtain 

dv^ -= / / 5 i^dvxdv2Xn (1^ 



7rV2^ Jr VTTv^ 7r^V2^ JrJr I + vf + v. 



These are of the order l/^/s as s — )■ 0. Now we indicate why the contributions 
from the poles of the A^^, when we shrink the contours, are of lower order. 

Consider the permutations (±2 1), when A^^ = S{^2,^i)- (Times ^f^ when a = 
(—2 1); this has no effect on what follows.) With the ^2-integration over Cr, we 
shrink the ^i-contour to Ci. Then when we shrink the ^2-contour we pass the pole at 
^2 = 2 — ^If ^ for all ^1 G Ci. The residue is a constant times 



Ci 






(19) 



where we replaced ^i by ^. With either branch of ^^/^ we may write 

(^_1)2 (^1/2 _ ^-1/2)2 



-1 



2^-1 2 - ^ 

On Ci the numerator is negative real (except when ^ = 1) while the denominator 
lies in the right half-plane. Thus the quotient lies in the left half-plane (except 
when ^ = 1) so we may take any s > in the integral. Since again the main 
contribution comes from near ^ = 1, we set ^ = e^^ and make the replacements 

2 — ^~^ -^ 2 — {1 — iv) = 1 + iv ^ e™, and we get 

^ ;i^(^-2/) dv = J~s I -^^ e^^^<^-y) dv. 



S + v'^ 7k 1 + f 2 

After the scaling this becomes independent of s. 

The factor {^ — 1)^ in the integrand in flT9|) was important. It also appears in 
the residues for the other integrals, which also become independent of s by similar 
computations. We omit the details. 

Thus when we scale the contributions from the poles of the A^ are independent 
of s, and so of lower order than the main terms (ITS]) . 

Set 

^iv{x-y) _j_ ^iv{x+y) 



Ji{x,y) = 


271 


M^^y) = 


1 

7r2 



vT+v^ 




;2 

gi {vi x~V2 y) 
2 



dv 



l + vf + v. 



dvi (if 2. 



^•^The integrals are Kq Bessel functions, but this fact is not useful. 
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We showed that l + 'jLl(s), when scaled, is equal to 2'j/\/2s {J1 + J2) plus an operator 
independent of s. Therefore we presume that to a first approximation the scaled 
operator (1 +'yL'l{s))~^ is equal to (27)""'^ v^ ( Ji + .12)^^- Also, we see from (TTOj) and 
the symmetry of p(x, y; s) that 

po{x) = p{x,0;s) r^ \ -e'"^"" as s -)■ 0. 
This gives the Conjecture 

((/ + 7i:?(s))-'Po,l)~^=((^l + ^2)-'e-^l) ass^O. (20) 

7V2s 

We shall show that the inner product equals 2/7r. Assume this, and the conjecture, 
for now. Then from ( IT6l) and the asymptotics p(0, 0; s) ~ a/2/s as s — > we find 
that 

E/' — ~~~~\ 1 <^^ -2 

X\<X2 



Therefore by the Tauberian theorem. 



1«2 



{VxiVxi) ^ 5-^ as t -^ 00. 

TT ^ 

Xl<X2 

Since 

{^N{tf) = J2 (Vx, Vx2) = 2 J2 (^-1 Vx^) + (AAr(t)), 

a::i,X2>0 a;i<a::2 

and AN{t) = o(t), we get the result stated in Section II, 

{AN{tf) r^ -^t as t ^ 00. 

n 7^ 

Now we show that the inner product in fl20l) equals 2/7r. We first use an observa- 
tion about Wiener- Hopf plus Hankel operators, of which Ji is one. (Recall that our 
operators act on functions on M+.) Suppose we have such an operator with kernel 
J{x — y) + J{x + y), so that the result of its action on a function / on M+ is 

{Jf){x) = / {J{x -y) + J{x + y)) f{y) dy. 
Jo 

If we extend / to an even function on M, then this equals 

J{x - y) f{y) dy. 
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So the operator becomes simply convolution by J . 

Also, it follows from the fact that J2(x, y) is even in y that if /, defined on R"*", is 
extended to be even on M then 

1 /""^ 

{J2f){x) = -j J2{x,y)f{y)dy. 



To recapitulate, define 



Ki{x,y) = — / -j===dv, 

^TT Jk VI + v^ 



i^2(a;,y) = T^-j / / ^ , ^ 2 I 2 '^^1 '^^2 



2^ ia iiR 1 + t;f + t;2 

acting on functions on M. Then on M+ we have (Ji + J2)/ = (-^1 + -^2)/; where the 
/ on the right is the even extension of the / on the left. 

It follows from this that for an even function / on M, ( Ji + 1/2)^^/ is the restriction 
to M+ of [Ki + -^'2)^^/0 Therefore the inner product in f l20|) . which is over M^, is 
equal to 

\{{K, + K,)-'e-\^\A), (21) 

where this inner product is over M. 

Because of the forms of the kernels of Ki and 7^2, the operators simplify when 
we conjugate with the Fourier transform. The operator Ki becomes i^i, which is 
multiplication by the function l/yl + i'^, and the operator K2 becomes 7^2, which 
has kernel 

K2{u,v) = -— — — — -. 
n 1 + u^ + V' 

Since the Fourier transform of 1 with factor l/27r outside the integral is 60 and the 
Fourier transform of e~'^' with factor 1 is 2/(1 + f ^), we see that fl2T|) (which equals 
the inner product in fl20l) ) equals 

{{K, + K2r\l + vY\6o). 
If we use 

{k^ + K^y = k-^/^ (1 + k-^'^ £2 k-^'\^ k-^'^ 

^ — 1/2 

and the fact that K^ 5q = 5q, we see that the above equals 

((/ + K)-V,^o), (22) 



^•^This uses that Ki + K2 commutes with the operator f{x) -^ /(—a;). 
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where 

If we conjugate with the unitary operator f{u) — )■ (coshx)^/^/(sinhx) we find 
that ([22]) equals ((/ + K)'^ ip, 6o), where 

~ , , , , \i/2 f>/ • 1 . , \ / , Ni/2 1 coshx coshy 
K{x,y) = (coshx) ' A (smhx, smhy) (coshy) ' 



From 



one sees that 



TT 1 + sinh^ X + sinh^ y ' 
tp{x) = (coshx)^'^'?/^(sinhx) = sechx. 

coshx coshy = -(cosh(x + y) + cosh(x — y)), 
1 + sinh^ X + sinh^ y = cosh(x + y) cosh(x — y) 



K{x, y) = -— (sech {x — y) + sech (x + y)). 

ZTT 



Now ip is even, and when K is restricted to the space of even functions it equals the 
operator with convolution kernel sech {x — y)/7f. Conjugating with the Fourier trans- 
form, this operator becomes multiplication by sech (7r,^/2). The Fourier transform 
of sechx with factor l/27r is (1/2) sech {tt^/2), and the Fourier transform of Sq with 
factor 1 is 1. Therefore 

2 J_^ 1 + sech (774/2) TT 

Thus the inner product in ( 120|) is 2/7r as claimed. 

IV. Appendix — The formulas for p(x, y;£) 

The Weyl group B„ is the group of signed permutations, functions a : [1, n] — )■ 
[—n, —1] U [1, n] such that \a\ G §„. An inversion in B„ is a pair (±cr(i), a(j)) with 
i < j such that ±o"(«) > o"(j). We write r = p/q. 

We define 

P + g^f-f 1-r^ 1 

and then define 

^a = JJ r{^^^i)) X JJ{S'(4'a, ^b) ■ {a,b) is an inversion in B„}, 

o-(j)<0 
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with the convention ^_a = t/^q. 

The formula, vahd when g 7^ 0, is 






(23) 



where x = {xi, . . . , x„} and y = {yi, . . . , x„}. The domain of integration is 

[J C/,(i) X ■ ■■ X C^(„), 

where the Ca are circles with center l/2g and distinct radii Ra- The i?^ should be so 
large that S{^,^') is analytic for ^, ^' on and outside CaL3 



In the special case n = 1, 



2vrz Jcfl 



(^•a;— «— 1 I > ^x /■—x—y—2 



MOt 



d^, 



where Cr is a circle around of radius R > 1. 

The formulas do not hold for p = 1 TASEP on Z+. But then the probability is 
the same as for TASEP on Z, and we have then [3] (or |5l p. 820]) 

p(x, y; t) = det fj (1 - ^y-^^^^-y.-^e'^iO dA , (24) 

where Cr is a circle with center and radius r < 1. 



•T 



Acknowledgments 

We thank Lauren Williams for elaborating for us some of the results of |1]. 

This work was supported by the National Science Foundation through grants 
DMS-1207995 (first author) and DMS-0854934 (second author). 



References 

[1] B. Derrida, M. Evans, V. Hakim, and V. Pasquier, Exact solution of a ID 
asymmetric exclusion model using a matrix formulation, J. Phys. A: Math. 
Gen. 26 (1993), 1493-1517. 



^^We cannot simply take C x ■ ■ ■ x C with C a circle with large radius, because then there would 
be nonintegrable singularities of S{£,a, ^b) on the contour when a > 0, 6 < 0. However by taking the 
Ra — >■ R we can interpret each integral as a symmetric distribution supported on C x ■ ■ ■ x C applied 
to the product in the integrand. 

21 



[2] P. L. Krapivsky, Symmetric exclusion process with a localized source, Phys. Rev. 
E86 (2012), 041103. 

[3] G. M. Schiitz, Exact solution of the master equation for the asymmetric exclu- 
sion process, J. Stat. Phys. 88 (1997), 427-445. 

[4] T. Sasamoto and L. Williams, Combinatorics of the asymmetric exclusion pro- 
cess on a semi-infinite lattice, larXiv:1204. 11141 

[5] C A. Tracy and H. Widom, Integral formulas for the asymmetric simple exclu- 
sion process, Comm. Math. Phys. 279 (2008), 815-844. 

[6] C A. Tracy and H. Widom, Erratum to "Integral Formulas for the Asymmetric 
Simple Exclusion Process", Comm. Math. Phys., 304 (2011) 875-878. 

[7] C A. Tracy and H. Widom, The Bose gas and asymmetric simple exclusion 
process on the half-line, J. Stat. Phys., 150 (2013), 1-12. 



22 



